Exploring Two Higgs Doublet Models Through Higgs Production 

Chien-Yi Chen and S. Dawson 
Department of Physics, Brookhaven National Laboratory, Upton, New York, 11973 



(Dated: March 4, 2013) 

Abstract 

We discuss the connections between the recently observed Higgs-like particle and rare B decays 
in the context of two Higgs doublet models (2HDMs). The measured decays of the Higgs boson to 
fermions and gauge bosons, along with the observation of the decay B s — > fi + fi~, place stringent 
restrictions on the allowed parameter space of 2 Higgs doublet models. Future measurements of 
— > 77 can potentially exclude type I 2HDMs, while the parameters of other 2HDMs are already 
severely restricted. The recent observations of the h° — > t + t~ and h° — > bb decays further constrain 
the models. 
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I. INTRODUCTION 



The discovery of a Higgs-like particle at the LHC with a mass m^o ~ 125 GeV is the 
beginning of the exploration of the source of electroweak symmetry breaking. The new 
particle has been observed in a number of different final states with rates roughly consistent 
with the Standard Model predictions, although the decay rate to 77 is slightly high. The 
crucial task now is to extract the properties (spin/parity and couplings to fermions and gauge 
bosons) of this particle as precisely as possible and determine whether they correspond 
to those predicted by the minimal Standard Model. It is also important to explore the 
possibility that there is a spectrum of Higgs-like states. 

Many beyond the Standard Model (BSM ) scenarios have a Higgs-like particle which has 
Standard Model -like couplings to fermions and gauge bosons in the low energy limit, and 
precision measurements of the Higgs particle couplings can serve to limit the parameters of 
BSM models. A well motivated extension of the Standard Model is obtained by adding a 
second SU(2) L Higgs doublet, leading to 5 physical Higgs particles: h°, H°, A , and H ± . 
We will consider the possibility that the particle observed at the LHC is the lightest neutral 
Higgs particle, h°, of a 2 Higgs doublet model (2HDM.) 1 

The 2HDM models generically have tree level flavor changing neutral currents (FCNCs) 
from Higgs exchanges unless there is a global or discrete symmetry which forbids such 
interactions and therefore we consider only the class of models where there is a discrete 



Z2 symmetry such that one of the fermions couples only to a single Higgs doublet. FCNCs 
are highly suppressed in this case and provide stringent limits on the parameters of the 
models. There are four possibilities for 2HDMs of this type which are typically called the 
Type I, Type II, Lepton Specific, and Flipped modelsj^]. We review the limits from FCNCs 
on these models, and determine the range of parameters consistent with the latest Higgs 
measurements at the LHC. Our goal is to study the extent to which LHC measurements in 
the Higgs sector can restrict the possibilities for 2HDMs. Previous works have examined the 
possibility of enhancing the branching ratio h° — > 77(1, 6, 7|, and various other channels js- 



12J in the 2HDMs and we discuss the implications of a measurement of h° — > 77 which 



differs from the Standard Model prediction and demonstrate the impact of the h° — >■ bb and 



1 The possibility that the observed particle is the heavier neutral Higgs particle, H Q , of a 2HDM has been 
examined elsewhere 111- 
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TABLE I: Neutral Higgs Couplings in the 2HDMs 
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TABLE II: Charged Higgs Couplings in the 2HDMs 





I 


II 


Lepton Specific 


Flipped 




cot/3 


cot j3 


cot/3 


cot/3 


Abb 


cot/3 


— tan j3 


cot/3 


— tan /3 


A-TT 


cot/3 


— tan /3 


— tan (5 


cot/3 



/z° — >• t + t measurements. Finally, we show how the recent measurement of B s — > fi + fi 



131 ] serves to further restrict 2HDMs. 



II. REVIEW OF 2HDMS 



The Higgs sector of the 2HDMs is parameterized in terms of tan/3 = where V\ and f 2 
are the vacuum expectation values of the two Higgs doublets and satisfy = g 2 {v f+v%)/2, 
and the angle a which diagonalizes the neutral Higgs mass matrix. The couplings of the 
lightest Higgs boson, h°, to the Standard Model particles are parameterized as: 

£ = -Z t g iih ^Jj t h - T Jt=w ^ z g hVV 2 -^V^h° , (1) 

v v 

where guh = ghvv = 1 in the Standard Model and v = 246GeV. The h° coupling to gauge 
bosons is the same for all four models considered here, while the couplings to fermions 
differentiates between the models. The Higgs Yukawa couplings normalized to their Standard 
Model values are summarized in Table [B The Standard Model couplings are obtained for 
sin(/3 — a) = 1, sin a = — cos (3, and cos a = sin/3. The charged Higgs-fermion couplings can 
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be written as, 
£ = 



V2M ] 



YV 



t \ tt m t P L - X bb m b P R bH 
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XnmiP R lH + + h.c. 



(2) 



ir 



where Pl,r = m i is the charged lepton mass and the coefficients A// are given in Table 

o 

Perturbative unitarity of the Yukawa couplings is violated if 

2 
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(3) 



For all models considered here this requires tan/3 > 0.28. Model II and the Flipped model 
require tan/3 < 140, while the Lepton-specific model has the limit tan/3 < 350 from per- 
turbative unitarity. Other bounds on perturbative unitarity have been derived by requiring 
that the quartic couplings in the scalar potential remain positive up to a high scale 14J and 
by considering the perturbative unitarity of gauge boson scattering 15j. These bound typi- 
cally give much lower upper bounds on tan /3 than those derived from perturbative unitarity 
of the Yukawa couplings, but we will not consider them further since these limits can be 
evaded by postulating new physics at some high scale. 

There are also strong bounds on the 2HDMs from precision electroweak measure- 
ments, which can be parameterized by the oblique parameters, S,T and U. In the limit, 



M H o, Ma, M H ± » M z and subtracting the Standard Model contribution 
1 
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For M A = M H o = M H ± and M A » M z , 

aAT ~ 

AS ~ -^-cos 2 (/3-a) 
YLtx 





5" 




" 6. 



(6) 



For = M#o = M#± ~ 1 TeV, AT provides no useful limit on sin(/3 — a), while AS 
requires cos(/3 — a) ~ 1. The most significant restrictions come from B decays, as discussed 
in the next section. 

The branching ratios to fermions in the 2HDMs are simply scaled from the couplings of 
Table |T] and the total h° decay width. We use the Standard Model branching ratios from 
the LHC Higgs cross section working group for M h o = 125 GeV j 
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(7) 



The total widths in the 2HDMs are given by, : 



Y"2HDM 

~ .577^ + .029^ cc + -063^ rr + .241sin 2 (/3-a) 

1 ho 



v(h° -> 77 ) 2 ™^ r(fr° -> ^) 2HDA/ 

+ r(/i°^ 77 ) SA/ + T(h°^gg) SM + "' () 

where we neglect other contributions which are smaller than the h° — > 77 branching ratio 
and assume that there is no new physics beyond the 2HDM. 

n □ 

The decay width to 77 is found using the exact form factors of Refs. [19|, |20J. For 



We assume guww = 9hzz- 
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Myp = 125 GeV, the dominant contributions are 



,0 _^ nrn M2HDM ( 2 

.2%g hu - .004g hbb - .003Qg hTT - 1.27 g hWW 



-> 77)] S 
[F(h° jj)] sm 



+{. 0057^ 66 +.0033^ TT | . (9) 
The tri-linear h°H + H~ coupling is not enhanced by a mass factor, and we neglect its neg- 



ligible contribution 



2l|. 



{ 1.06^-. 06^ J +|.086^ 66 | . (10) 



Similarly, the decay to gluons proceeds predominantly through top and bottom quark 
loops, 

T (gg -»■ h 0fHDM 

T(gg h°)S" 

For large tan/3, the b quark loop can contribute significantly to the gg — > h° production 

channel in Model II and the Flipped model. In our numerical studies, we include both the 

b and t contributions exactly in all cases . 

Using the results given above, the production from gluon fusion and the following decay 

to 77 or to // pairs can be described by, 

F [a{gg -» h°)BR(h° -> 77 )] 2 ^ 
77 " [a(gg -> h°)BR(h° -> 77 )] 5M 

^ = W(gg -> h°)BR(h° -> /7)T DM fnl 

J/ [a(^ -> h°)BR(h° -> //)] SM ' 

The Higgs signal from gluon fusion with the subsequent decay to 77 is shown in Fig. [1] 
for v^S* = 8 TeV 3 . It is interesting to note that in Model I, it is not possible to obtain R 99 ^ 
larger than ~ 1.2, so a future measurement which confirms the current deviation of R 99 ^ 
from 1 could serve to exclude Model I 4 . For R 99 f > 1, only a narrow region of a ~ —0.5 
can be obtained. The Lepton Specific model is similar to Model I at small tan/3, while at 
large tan /3 ( for fixed non-zero a) the total width is enhanced by the large branching ratio 
to t + t~, which decreases R 99 ^ . A value of R 99 ^ ~ 1.2 is only consistent with the Lepton 
Specific model for a ~ and tan /3 > 8. Model II and the Flipped model have similar 
predictions for R 99 ^ with a wide range of a and tan/3 values consistent with R 99 ^ ~ 1. 

The gluon fusion production with decay to t + t - is shown in Fig. [2j In Model I and the 
Flipped model, the Standard Model rate can be obtained for small a, relatively indepen- 
dently of tan /3. Model II and the Lepton Specific model only approximate the Standard 



There is very little change going from = 7 TeV to VS = 8 TeV. 
4 This has also been noted in Ref. [31. 
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FIG. 1: The contours correspond to constant E^f as defined in the text for \/~S = 8 TeV. 
For Type I (a): The blue (dashed), black (solid), and red (dotted) lines are Wyf = .9,1,1.1, 
respectively. For Type II (b): The blue, black, and red lines are B3f F = .9, 1., 1.4, respectively. For 
the Lepton Specific model (c): The blue, black, and red lines are RHjHf = .9,1., 1.2, respectively. 
For the Flipped model(d): The blue, black, and red lines are Rlyy F = .9, 1., 1.4, respectively. 

Model rate for very specific values of a and tan /3 and hence a precise measurements of R 9 T 9 T F 
can serve to restrict these models significantly. 

Preliminary LHC measurements of h° — > bb and h° — > t + t~ in the vector boson fusion 
(VBF) and Vh° channels are particularly interesting and will serve to distinguish between 
the different 2HDMs. Both VBF and Vh° production follow the same simple scaling between 



7 



10 



50, 



1 
1 
1 
1 
1 
1 


1 


1 

j 0.9 

1 
1 
1 


1 
1 




i 1 
\ 1 
\ \ 
\ \ 






■\\ 




1.1 \\\ 

















-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



(a) 



(b) 



10 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



(c) 




(d) 



FIG. 2: The contours correspond to constant R%% F as denned in the text for y/S = 8 TeV. 
For Type I (a): The blue (dashed), black (solid), and red (dotted) lines are B%- 9 T F = .9,1,1.1, 
respectively. For Type II (b): The blue, black, and red lines are Rt% F = -9, 1., 1.4, respectively. For 



the Lepton-Specific model (c): The blue, black, and red lines are R 



ggF 



.9,1., 1.2, respectively. 



For the Flipped Model(d): The blue, black, and red lines are R% gF = .9, 1., 1.4, respectively. 

the Standard Model rates and the 2HDM predictions, 

sin 2 (P- a)g 2 hff . (12) 
From Fig. [31 it is clear that the Standard Model rates, R^ b BF,Vh ~ 1 only occurs for very 



r>VBF,Vh ■ 2 1 

R ff 



specific values of a and tan/?. In Model I and the Lepton Specific model, R^ b BF,Vh ~ 1.1 
requires tan < 1, while Model II and the Flipped model require a ~ —0.5 and tan /3 > 4 — 8 
for R^ b BF,Vh ~ 1.4. The rates to t + t~ normalized to the Standard Model predictions are 
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FIG. 3: The contours correspond to constant R^ BF ' Vh as denned in the text. For Type I (a): 
The blue (dashed), black (solid), and red (dotted) lines are R^ b BF,Vh = .9, 1, 1.1, respectively. For 
Type II (b): The blue, black, and red lines are R^ b BF,Vh = -9, 1., 1.4, respectively. For the Lepton- 
Specific model (c): The blue, black, and red lines are R^ b BF,Vh = .9, 1., 1.2, respectively. For the 
Flipped Model(d): The blue, black, and red lines are R^ b BF,Vh = -9, 1., 1.4, respectively. 



shown in Fig. gj Note that in Models I and II, R bb BF ' Vh = R v T BFyh . This is because the 
h° coupling to bb is the same as the h° coupling to t + t~ in these two models, as shown in 
Table [H As a result Fig. is identical to Fig. H^a) and similarly Fig. [3](b) and Fig.H^b) are 
also identical to each other. One can also see that Fig. 3(a) and Fig. 3(c) are very similar 
since the coupling ghbb for Model I and the Lepton-Specific model are the same. They are 
not identical because the total widths in these two models are different. This similarity also 
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(c) (d) 

FIG. 4: The contours correspond to constant Rr^ F ' Vh as denned in the text. For Type I (a): 
The blue (dashed), black (solid), and red (dotted) lines are Rr^ F,Vh = .9, 1, 1.1, respectively. For 
Type II (b): The blue, black, and red lines are Rr^ F,Vh = .9, 1., 1.4, respectively. For the Lepton 
Specific model (c): The blue, black, and red lines are R}f-^ F,Vh = .9, 1., 1.2, respectively. For the 
Flipped Model(d): The blue, black, and red lines are Rr^ F,Vh = .9, 1., 1.4, respectively. 

appears in Fig. 3(b) and Fig. 3(d) because the h° coupling to bb is the same for Model II 
and the Flipped model. 

III. REVIEW OF LIMITS FROM FLAVOR PHYSICS 



Limits on 2HDMs have been examined by many authors and we briefly update the results 
of Ref. 22j using the Superlso program 23M25| in order to examine the restrictions on 
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FIG. 5: Limits on 2HDMs from AAf^ d . Regions below the blue and red regions are excluded at 2 
and 3a, respectively. 

tan /3 in the 2HDMs which are relevant for Higgs coupling studies. Small values of tan /3 , 
tan/3 < 0.35, for Mh+ < 2 TeV are excluded at 3a in all 2HDMs considered here by the 
experimental measurement of AM^ 26| . 



AM 



B d \exp 



0.507 ±0.004 ps' 1 . 



(13) 



The limits from AMg d are identical in all 2HDMs studied in this work and the exclusion 
regions are shown in Fig. [5J For M H + ~ 300 GeV, tan/3 < 1.03 is excluded at 3a. (For 
Mh+ ~ 2 Te^, tan/3 < .35 is excluded at 3a.) 

Stringent bounds also arise from the experimental measurement of B — > X s ^ 



26 



27] 



BR(B X s ~?) 



exp~ 



(3.55 ±0.24 ±0.09) x 10" 



(14) 



and are shown in Fig. [6|22 



For M H+ = 1 TeV, tan/3 < 0.69(0.16) is excluded 
at 3a in Type I and the Lepton Specific model (Type II and the Flipped Model). For 
M H + = 300 GeV the restrictions are still fairly stringent: tan/3 < 1.36(0.88) is excluded at 
3a in Type I and the Lepton Specific model (Type II and the Flipped Model). For all values 
of tan/3, M H + < 300 GeV is excluded in the Type II and Flipped models. The numerical 
results agree with Refs. 



22 



29 



31| in the appropriate limits. 



Finally, the recent observation of the decay B s — > 13[ has important implications 
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FIG. 6: Limits from the measurement of B — > X s j in Type I and the Lepton Specific model (a). 
The regions below the blue (upper) and red (lower) curves are excluded at 2 and 3a. The limits 
in Type II and the Flipped Model are in (b) and the regions to the left of the blue and red curves 
are excluded at 2 and 3a, respectively. 



for the 2HDMs, 



BR(B 8 ->■ 



I exp 



(3.2+H) x 10' 



in good agreement with the Standard Model prediction 32| 



BR(B S -»■ fj+fi 



\SM Z 



(3.23 ±0.27) x 10" 



(15) 



(16) 



New physics effects in B s — > fi + fi come predominantly from the charged Higgs exchanges. 
The contributions to B s — > in the type II two Higgs doublet model have been computed 



in the large tan j3 regime in Refs. 



33 



3M35| . and adapted to the general Two Higgs doublet 
models considered here in the Superlso program 24J. 5 

The branching ratio for B s — > fi + fi~ in the Type I model is shown in FigJT^a) as a function 
of Mh+ for various values of tan/3. It is apparent that for Mh+ > 500 GeV, the branching 
ratio is almost a constant, independent of tan/3. In FigJT^b), the excluded region is shown 



It is possible that there are contributions which are not enhanced by large tan/3 in the Type II model 
which could be relevant in the Type I, Lepton Specific or Flipped models. Such contributions are not 
included in the Superlso program. 
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(a) (b) 

FIG. 7: Branching ratios for B s -> fj+fi- in Model I (a) with M H o = M A = 300 GeV for tan /3 =5 
(blue), 15 (red), and 30 (green). The excluded region from B s — > in Model I is shown in 

(b) with M H o = Ma- The red (blue+red) region of (b) is excluded at 3(2)<r for a = —0.88 and 
tan (3 = 1. These values of a and tan j3 correspond to the best fit to the Higgs data which is derived 
in Section IV. This plot, however, is rather insensitive to the precise value of tan f3. 

for the parameters which best fit the Higgs data (derived in the next section). However, for 
M H + and Mh > 500 GeV, the excluded region is not sensitive to the value of tan/3. Even 
for smaller values of Mh+ and Mjj, the sensitivity to tan/3 is small. 

In the Type II model, the branching ratio B s — > has a significant dependence on 

tan /3 for small Mh+ , and goes to a constant for very large M#+ . In this model high values 
of tan /3 are excluded at 3a for small M^+, while there is a 2a excluded region at large tan {3, 
as shown in Fig. [SJ The dependence of the excluded region on the choice of neutral Higgs 
masses, M^o and Ma, is shown in Fig. [9] for a = —0.02 and tan/3 = 60, (which corresponds 
to the best fit to the Higgs data derived in Section IV). 

For the Lepton Specific case, the dominant contribution to B s — > is proportional 

to AftA^, and so the branching ratio is insensitive to tan/3. In Fig. [TO], we show the 
BR(B S — > /U + /i~) as a function of and the regions which are excluded at 2 and 3a from 
this decay. We see that the excluded region does not depend on tan /3. For heavy charged 
Higgs masses, the branching ratio approaches a constant. The dependence on the choice of 
neutral Higgs masses is shown in Fig. (TTJ From Fig. [TOT a), it is clear that the branching 
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FIG. 8: Branching ratios for B s -> in Model II (a) with M H o = M A = 300 GeV for 

tan/? =5 (blue), 15(red), 30 (green), and 60 (black). (At high M H +, tan/3 increases going from 
the top to bottom curves.) The excluded region from B s — > [J> + fJ>~ in Model II is shown in (b) with 
M H o = M A = 145 GeV. The red (blue+red) region of (b) is excluded at 3(2)a. 

ratio for B s — > increases for M H + > 150 GeV as M H o and Ma decrease. This leads to 

the exclusion of the region between M H o = Ma < 100 GeV and M H + > 150 GeV at 2 a in 
Fig. dU 

The branching ratio for B s fi + fi~ for the Flipped model is shown in Fig. [T2] for several 
values of tan/3 and for M H o = Ma = 300 GeV. For M H + > 300 GeV, the branching ratio 
is insensitive to the input parameters. 



IV. RESULTS FROM HIGGS MEASUREMENTS 

We do a simple x 2 fit to the data shown in Tables HTT1 and HVl assuming M^o = 125 GeV. 
We follow the standard definition of \ = ' ^^asy — — , where R IHUM represents predic- 
tions for the signal strength from the 2HDMs and R mcas stands for the most recent results of 
the measured signal strength shown in Tables HTT1 and HVl by the ATLAS and CMS collabora- 
tions at the LHC. <j meas denotes the uncertainty of R meas . 6 The results (with no constraints 

6 The VBF tagged channels have a small contribution from gluon fusion. The ATLAS results (Table 1 
of Ref. [i^) explicitly separate the true VBF contribution from the gluon fusion channel using Monte 
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FIG. 9: Excluded region from B s — > fi + /i~ in Model II with M H o = Ma for a = —0.02 and 
tan/? = 60. (These values of a and (3 correspond to the best fit to the Higgs data which is derived 
in Section IV.) The red (blue+red) region of is excluded at 3(2)a. 

from flavor physics) are shown in Fig. O In all cases, the \ 2 minima occurs for a ~ 1.3 — 1.4 
and tan/3 ~ 0.2 — 0.3. The results of Sec. 1111} however, show that even for heavy Mjj+ and 
heavy M#o and Ma, such small values of tan /3 are not allowed in the 2HDMs we consider 
here due to constraints from the B sector. 

We perform a constrained fit to the data requiring tan /3 > 1, which is consistent with B 
physics data of the previous section. In all models, the results of Fig. [T3] show that there 
are large regions of parameter space allowed at both the 2 and 3a confidence levels and the 
results have only a mild dependence on tan (3. In Model II and the Flipped model, a = is 
not allowed, primarily due to the h° — > bb measurement. In Model II, the Lepton Specific, 
and the Flipped model, only a small range of a is allowed and the \ 2 minimum occurs for 
large tan/3 (60,54, and 77, respectively). In Model I, a fairly large range of a is consistent 
with the data and the x 2 minimum is at tan /3 = 1. 

Carlo. The CMS VBF results of Ref. [3] (Table 1) contain a 30 — 50% contamination from gluon fusion, 
estimated from Monte Carlo. We assume a 30% contamination of the CMS VBF result from gluon fusion, 
although the results of Fig. [13] are not sensitive to this assumption. 
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FIG. 10: Branching ratios for B s — > in the Lepton Specific model (a). (The branching ratio 

is almost independent of tan/3.) The excluded region for Mjjo = Ma = 145 GeV is shown in (b). 
The red (blue+red) region of is excluded at 3(2)er. 

V. CONCLUSIONS 

We have considered four variations of 2HDMs which have a Z 2 symmetry suppressing 
tree level FCNCs. Higgs production and decay in the 2HDMs can be significantly different 
than in the Standard Model and only small regions of a — tan /3 can produce rates which 
are consistent with the experimental results from the LHC. Further, the parameters of these 
models are strongly constrained by measurements in the B sector. In particular, limits on 
AM# d require tan /3 > 0.35 for Mh+ < 2 TeV in all 2HDMs considered here. For each model, 
we have also shown the regions in parameter space which are allowed by the measurement of 
B s — > fi + fi~ for the parameters which correspond to the best fit to the Higgs data. Unitarity 
also restricts the allowed regions to have tan/3 > 0.28. Our major result is shown in Fig. 
[TBI where we show the regions of a — tan /3 which are consistent with the Higgs cross section 
and branching ratio results at the 2 and 3a level. None of the models we studied can be 
excluded by current measurements. 
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FIG. 11: Excluded region from B s — > [i + [i~ in the Lepton Specific model with Mjjo = Ma 
and a = —0.02, the best fit to the Higgs data, derived in Section IV. (The excluded region is 
independent of tan/3.) The red (blue) region of is excluded at 3(2)<r. 
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